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The theoretical description of interacting fermions in one spatial dimension is simplified by 
the fact that the low energy spectrum of noninteracting fermions is identical to the one for a 
harmonic chain. This fermion-boson transmutation allows to describe interacting fermions as a 
system of coupled oscillators. Tomonaga's model of interacting nonrelativistic fermions on a ring is 
presented and first discussed in low order perturbation theory. After introducing the concept of two 
independent species of right and left moving fermions the exact solution of the Tomonaga-Luttinger 
■ model is discussed in detail. The momentum distribution and spectral functions are calculated using 

| the method of the bosonization of the field operator. The general Luttinger liquid phenomenology 

is shortly discussed. 
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I. INTRODUCTION 



The effects of the Coulomb interaction of electrons in metals can only be described approximately!!!. With the 
progress in producing artificial low dimensional structures the theoretical work on one- dimensional interacting fermions 
• ■ has gained importances. Special features of the spectrum of low energy excitations in one dimension allow exact 
solutions of models of interacting fermions. The main idea can already be understood by working with noninteracting 
fermions, which have the same spectrum of excitations as a harmonic chain.-. This "fcrmion-boson-transmutation" 
(FBT) was discussed in a recent publication on different levels of sophistication!! In the following this paper is referred 
*^ . to as I. The basic idea can be understood on the level of elementary quantum statistical mechanics. Only in chapter 
V. of I the method of second quantization was introduced to calculate the momentum distribution of noninteracting 
fermions in the canonical ensemble. It was pointed out that the technique presented can also be used in the description 
of interacting fermions. From the response we obtained to I it is clear that many readers would have liked to see more 
explicitly the application to this problem. The purpose of this paper is to fill this gap on a level which requires only 
i a basic knowledge of the method of second quantization but none of relativistic quantum field theory 

The exact solution to the "Tomonaga-Luttinger (JL-) model" is presented starting with Tomonaga's model of 
spinless interacting nonrelativistic fermions on a ringa. In section II the momentum distribution in the interacting 
' ground state is calculated in leading order perturbation theory. This already demonstrates one of the important ideas 
towards the exact solution and on the quantitative level shows a sign of the infrared singularities later found in the 
t— I ' exact solution. Tomonaga's concept of "right" and "left" moving electrons later extended by Luttingera is introduced 
, and the commutation relations of the Fourier components of the operator of the particle density arc derived. In order 
to avoid mathematical subtleties we work with two cut-offs - one for the number of momentum states ( "band cut-off" ) 
and one for the inverse range of the two-body interaction ( "interaction cut-off" ) . This allows to smoothly go from 
Tomanaga's original model to the model discussed by Luttinger. In section III. the exact low energy eigenvalues of the 
TL-model are obtained using the bosonization technique. As the most important "Luttinger-liquid" (LL) feature the 
momentum distribution in the ground state is calculated. Its derivative at the Fermi points diverges in a power law 
fashion determined by the important constant called "anomalous dimension" . One effort to experimentally verify LL- 
properties is photoemission from quasi-one dimensional conductorscJ . As the theoretical description of photoemission 
O ■ involves spectral functions they are discussed in section IV. The model including spin is introduced in section V. After 
the definition of "charge" and "spin" bosons the Hamiltonian of the model can be written as a sum of two commuting 
terms of the spinless type leading to the phenomenon of "spin-charge-separation". In section VI. the general LL- 
concept is presented and it is shortly discussed how the LL parameters which completely determine the low energy 
physics also for more complicated Id lattice models can be extracted from finite size numerical data. In appendix A 
the proof of the Kronig-identity which provides the bosonization of the kinetic energy is given. The bosonization of 
the field operator is discussed in appendix B. A novel straightforward approach is presented to determine the particle 
number changing contribution. We put h = 1 in this paper. 



II. INTERACTING FERMIONS: THE TOMONAGA MODEL 



The method of bosonization which was presented in I in the context of noninteracting electrons in one dimension is 
the key concept to understand ground state properties and the spectrum of excited states with low excitation energy 
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also for interacting fermions. In a seminal paper TomonagaS studied nonrelativistic fermions in the high density limit, 
in which the spatial range of the two-body interaction is much larger than the interparticle distance. He studied 
spinless fermions on a ring, i.e. using periodic boundary conditions. The Hamiltonian H is a sum of the kinetic energy 
T = yij—i pj /2m and a two-body interaction V 

v- = (i) 
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i /-V 2 /- L / 2 

- / da: / dx'^(x)^(x')V{x-x')il){x')il){x), 



L/2 J-L/2 



where in the second line we have switched to the language of second quantization, with tp(x) the field operator which 
annihilates a fermion at thfi-jpasition x. It can be expressed in terms of the annihilation operators c n of momentum 
states in the standard wayBEE3 

<P(x) = -±=J2 eiknXc "> ( 2 ) 

* n 

where the sum runs over momenta k n = 2nn/L 1 n 6 Z compatible with the periodic boundary condition. The ^(cj^) 
obey canonical anticommutation relations 

[Cm, Cn]+ = , [c m , C,j,] + = 5 mn L (3) 

In order to be consistent with the periodic boundary condition ip(x) = tj)(x + L) it is useful to expand the two body 
potential V(x — y) in a Fourier series 

V(x-y) = jY, etkniX ~ V)v ( k ^> W 

n 

r L/2 

v(k n ) = / V(x)e- iknX dx. 

J-L/2 

Special choices of the Fourier transformed potential v(k) will be discussed later. In case V^O) exists, we can express 
V in terms of the operator of the particle density 

N 

P( x ) = ^2 S ( X ~ Vi) = ^(x)ijj(x), (5) 
i=l 

where the Xi are the position operators of the particles. Using [ip(x),ip'(y)]+ — S(x — y) for x,y £ {—L/2, L/2] one 
obtains 



1 rL/2 ,L/2 i 

V = - I dx dx'V(x - x')p(x)p(x') - -V(0)Af, (6) 

-L/2 J-L/2 2 
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where we have introduced the particle number operator Af = J_ L / 2 p( x )dx. As in Tomonaga's paper we present the 
main results for the spinless model and discuss the modifications due to the electron spin later. It is useful to Fourier 
decompose the density operators 

/S(z) = ~£/3„e <fc »*, (7) 

n 

,L/2 

pn = / p{x)e~ lknX dx. 

J-L/2 

If we use the second quantized version p(x) = ip> (x)ip(x) the p n can be expressed as follows 

Pn =X] C n' c n'+n- ( 8 ) 
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In terms of the p n the interaction term is given by 

V = j ^2 v(k n )p n p^ n 



n>0 



-j^(0) \v{0)Af. 



(9) 



We have used v(k n ) = v(—k n ) as V(x) is real and [p n , p m ] = 0, which follows most easily from the fact that only the 
commuting position operators Xi enter in the definition of p(x) in Eq. (5). 

Now following Tomonaga we assume that the spatial range of the two-body interaction is much larger than the 
mean particle separation which is inversely proportional to the Fermi momentum hp. Two examples for v(k) with a 
finite value of t>(0) which will be used later in the paper are shown in Fig. 1. 




FIG. 1. Two forms of the dimensionless Fourier transform F(k) = v(k) / '(tvvf) of the two body potential. 

The ground state of the noninteracting system with TV fermions (N odd) is the Slater determinant \F > ("Fermi 
sea") in which the momentum states from — hp — —2nnp/L to kp are occupied, where np = (N — l)/2 and kp is 
the Fermi momentum. Ajiirst insight into the nature of the interacting ground state can be obtained using first order 
perturbation theory in Vu 



l^o >= \F> 



Q 



E 



(0) 



-V\F > +. 



T 



(10) 



where Q = 1 — \F >< F\ and E^ is the kinetic energy of the Fermi sea. Due to the projector Q only the first term 
on the right hand side (rhs) of Eq. (9) contributes to the first order correction in the above equation. As we assume 
v(k n ) to be zero for k n (much) larger than a cut-off k c <C kp only states p n p- n \F > with n of the order n c — k c L/2ir 
or smaller occur in V\F >. The second quantized form of p n Eq. (8) shows that p n p- n \F > is a two particle-two 
hole state. For n c > n > the density operator p_„ creates a particle-hole pair around the right Fermi point k = kp, 
while p n creates a particle- hole pair with opposite momentum around the left Fermi point k = —kp 



QV\F >= j v(k n 

n>0 



E E c t 

' — — rip- —n n'— + 1 



(ii) 



The kinetic energy of the two-particle-two hole state c n ,, c n »+ n c J n , c n >- n \F > is given by E^"' +e n » —£„"+„+£„' — e n 
where e n = (2irn/ L) 2 /2m e . As the particles and the holes are in the vicinity of the Fermi points one can linearize 
the energy dispersion around the two Fermi points in order to calculate e n > — t n '_ n and e n " — t n "^- n . With the Fermi 
energy ep = k F /2m e and the Fermi velocity vf — kp/m e one uses instead of the one parabolic dispersion two linear 
branches 



.(o) 



e„ !+ = ep + vp(n — np)2n /L n > 
e„._ = ep — vf(u + np)2ir/L n < 



(12) 
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As will be discussed later the use of the linear dispersions, also shown in Fig. 2, leads to an error which vanishes 
in the high density limit k F /k c — > oo. It simplifies the expression for the kinetic energy of the two particle-two hole 

state to + 2k n VF and leads to 

1 (h \ — n F +n 

| 0Q>= | F> _1^|M £ £ c t„ Cn „ +nC t, c „,_„|F>+... (13) 

n>0 n" ——np—n n'-n^+l 

Higher order perturbation theory yields additional terms with multiple particle-hole pairs. 

An important way to characterize the interacting ground state is to calculate its momentum distribution 
< 4>o\n m \(f>() >=< 4 > o\c T n c m\ l fio >, i-e. to study how the Fermi step distribution of the noninteracting electrons is 
modified by the interaction. It follows from Eq. (13) that the leading correction to the Fermi function is quadratic in 
the two-body potential. For m > n F one obtains e.g. 

<*.*»- h £ if (M) 2 < (14> 

n>m—riF n" — ~np—n 

where we have used that the two particle- two hole states are eigenstates of n m . The sum over n" can be trivially 
performed and yields a factor n. The sum over n is especially simple for the step potential v{k) = u(0)8 (fc 2 — fc 2 ) 
where O(-) is the unit step function. If the sum over 1/n is approximated by an integral one finally obtains 



k m > k F : <n m > 



<*>= - Q^-) 2 In (^^) ■ 6 [ kc - (k m - *,)] . (15) 



The logarithmic singularity for (k m — k F )/k c — > is an artefact of the perturbation expansion, as the occupancies have 
to be smaller than one. The exact solution for the momentum distribution < n rn > will be derived in detail later. In 
the thermodynamic limit L — » oo, N/L = const one obtains a continuous function n(k) which in the neighbourhood 
of k F is given by 



k — kr. 



sign (kp — k) . (16) 



The constant a is called the anomalous dimension. For weak two-body potentials a is small and one can use 
x a = exp(alna;) = 1 + a In a; + 0(a 2 ) for fixed x. By comparison with the perturbational result Eq. (14) one finds 
c = 1/2 and 



a 



(2) 



um 2 (it) 



2 \2ttv f J 



as the lowest order perturbational result for the anomalous dimension. A quantitative comparison of < h m >( 2 ) with 
the exact result will be given later. 



Perturbation theory shows that there are three types of electrons: Those with momenta k w kp, which move 
with velocities w vf which will be called "right movers", those with momenta k rj — k F and velocities « —v F ("left 
movers") and the "inert" electrons deep inj-the Fermi sea, which play no role in the "low energy physics". This basic 
insight of Tomonaga was used by Luttingem to further modify the model. In addition to the linearization Eq. (12) he 
extended the range of the possible k- values for both the right and left movers to 2irn/L with n positive and negative 
integers. This further simplifies the dispersion of the electrons, but introduces infinitely many negative energy electron 
states, which have to be filled in the ground state. As this infinite "Dirac sea" leads to various mathematical subtleties 
we prefer to work with a finite band cut-off. We introduce the momentum ks = 2-kuio/L > and allow states with 
momenta > —kg for right movers and < kg for left movers. In Fig. 2 we have chosen kg — 1.5k F and show the states 
added as the dot-dashed lines. 
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FIG. 2. Energy dispersion as a function of momentum. The dashed curve shows the usual "nonrelativistic" dispersion used 
in the original Tomonaga model and the full curve the linearized version Eq. (12). The dot-dashed parts are the additional 
states for ks = 1.5&F. The model discussed by Luttinger corresponds to ks — > oo. 



Later we also consider the limit ks — > oo. The kinetic energy of the model solved exactly in the next section reads 

T T L = v F k n4i, + c n,+ + ^ (~ V F k n) C n- C n - , (18) 

riG/ + n€i_ 

where I + = {—mo, —mo + 1, ...oo} and J_ = {— oo, ...mo — l,mo}. The creation and annihilation operators for the 
right (a = +) and left (a = — ) movers obey the fermion anticommutation relations (m £ I a ,n € Ip) 



-m,/3 



— <5n,m^a,/3l 



In the interaction term V Eq. (9) we replace p n p-n by (p n ,+ + Pn,-) (p-n,+ + P-n,-) where 

Ea t 
w n',n'+n C n' ,a c n'+n,a, 

n> 

with 



(19) 



(20) 



u n' ,n'+n 



1 for n', n' + n E I a 
else 



are the "density operators" of the right and left movers. Because w" m is symmetric in n and m the relation 

Ptn,a = P~m,a holds. 

It is not more difficult to work with the following generalization 

= J Y {gi{K) (Pn,+ P-n,+ + Pn - P-n-) + 9l{k n ) (/5„ i+ p- n - + p n - p- n ,+)} 



L 

n>0 

+^ 34 (0) {Ml + ^ 52 (0) • 2W + W_ 

where the Af a are the particle number operators 

M a = J]] c ln,a C m,c 



(21) 



(22) 



and we have dropped the terms proportional to the particle number operators as they can be later included in the 
chemical potential . The second line in Eq. (21) can also be expressed in terms of the total particle number operator 
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Af = Af+ +Af- and the "current operator" J = Af+ — AfJB. For the interaction functions in Eq. (21) we have used 
the standard notations. If we put gi(k) = gi (k) = v(k) the Tomonaga-Luttinger Hamiltonian Htl = Ttl + Vtl 
yields in perturbation theory in Vtl the same result for < fi m , + > as in Eq. (15). In Q Vtl\Ftl > only the terms 
proportional to /?„._ f>- n .+ \FTL > contribute for n > 0. Concerning the low energy physics the original model with 
the quadratic dispersion is equivalent to the Tomonaga-Luttinger (TL) model described by the Hamiltonian Htl in 
the high density limit kp/k c — > oo. 

There is an essential difference in the commutation relations of the p m ,a and the p rn . While as discussed previously 
[pm,Pn] — 0, which can be seen most easily in the first quantized version, the p m ,a do not all commute. From the 
anticommutation relations for the fermions on different branches it follows that [p m ,+ , Pn.-] — 0. The commutators 
\fim,oti Pn,a] on the other hand are nontrivial 

[Pm,a, Pn,a] = [ C m' ,a c m'+m,a, C^, a C„'+„, Q ] W m , iTn i+ m W n , „/ +n (23) 

m' 7 n' 

= (fin',m'+m c m',a c n'+n,a ~ 3m',n'+n C^/ Q C m '+ mjQ ^ w m ',m'+m w n',n'+n 

m'n' 

^ ] ( w m' ,m'+m ~ w m' ,m'+n) w m' ,m'+m+n C m' ,a c m'+m+n- 
m' 

In going to the last equation the summation index was changed from n' to m! in the second term. It is easy to 
see that the rhs of Eq. (23) vanishes if m and n have the same sign. Take, for example, right movers and m 
and n positive. Then due to the factor to*, m > +m+n the rhs of Eq. (21) can only be nonzero for m' > —m and 
m! + m + n > —mo. But then also ml + m > —mo and ml + n > —mo and the commutator vanishes because of the 
factor iu+, m , , m — w^, m , +n = 0. Next we consider the case n = — m shortly discussed in section V. in I. For m > 
one obtains for the right movers 

-mo+m-l 

[p m ,+,P- m ,+] = 4n',+ °m',+ =A/"+(m). (24) 

m' = — 7710 

We see that the commutator is given by the particle number of the lowest m right moving one-particle states. As 
explained earlier an interaction v(k) which drops off rapidly for |fc| > k c does not produce holes deep in the Fermi sea 
in the interacting ground state and eigenstates with low excitation energy. In this subspace of the total Hilbert space 
it is an excellent approximation which becomes asymptotically exact for mo — > oo to replace A/+(m,) by ml on the 
rhs of Eq. (24). In the interaction term Vtl Eq. (20) only p n>a for |n| -C np contribute. Therefore it is sufficient to 
discuss the commutation relations Eq. (23) for m and n with different sign for |m|, \n\ <C np. For m > — n > one 
obtains for right movers 

— mo— Tl— 1 —ino— 71— 1 

[Pm, + j P—n,+] = ^ t c m',+ c m'+7n—n,+ = ~ ^ ' c m'+m-n,+ c ln', + ' (25) 

The operator c| n , , tries to create an electron deep in the Fermi-Dirac sea, which is not possible in the space of low 
lying excited states, as all states at the bottom of the band are occupied. Therefore the commutator vanishes in this 
subspace. Analogous arguments apply to the left movers and we finally obtain 

[Pm,Q 5 Pn, p] = am 1 6 (26) 

The importance of this relation was first realized by TomonagaS who worked "udth ks = 0. Luttingeri made an error 
with the commutation relation for fcs = oo later corrected by Mattis and LiebLl. 

III. EXACT SOLUTION OF THE TL-MODEL USING BOSONIZATION 

The commutation relations Eq. (26) for the operators p m , a apart from a proper normalization factor, look like 
boson commutation relations. If we define the operators 

K = i J Am- fi»n>o (27) 

VW\ I Pn <- n K 
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and the corresponding adjoint operators tf n using p^ n a = p- n ,a Eq. (26) implies the usual bosonic commutation 
relations 

[b n ,b m ]=0, [b n , tfj = S mn l. (28) 

The kinetic energy Ttl in Eq. (18) consists of two commuting terms, which both are of the form discussed for fixed 
boundary condition in I. Therefore they both can be expressed in terms of the boson operators with the help of the 
Kronig identity^ which was stated in I without giving an explicit proof. This is therefore presented in appendix A. 
It yields 

t TL = y, »f (y) h 6 « 6 « + «*■ (i) + const - N ( 29 ) 

As we dropped terms proportional to M in Vtl which only lead to a renormalization of the chemical potential we 
also drop the term proporlional to Af in the kinetic energy in the following. Neglecting also a constant term in the 
total Hamiltonian it readscj 

Htl = £ \k n U + *M) (blb n + bl n b. n ) + (ftt bl n + b_ n b n ) ) 

n>0 ^ ^ n ' 7T J 

+ 2T ^ VNAf2 + VjJ ^ = Hb + H *J (30) 



with the velocities 



g4(0)+g 2 (0) ,„ n 
v N =v F ^ (31) 



Vj = Vp 



2tt 

.94(0) - g2 (0) 
2tt 



Here ujv determines the energy change for adding particles while vj enters the energy change when the difference in 
the number of right and left movers is changed. For our original model with (74(g) = 32(9) = v(q) the velocity vj 
equals the Fermi velocity vf of the noninteracting system. As the particle number operators Af a commute with the 
boson operators b m {bl rn ) the two terms Hg and H^- » in the Hamiltonian commute and can be treated separately. 

The boson part Hg of the Ti-Hamiltonian is bilinear in the boson creation and annihilation operators. As shown 
below it can be brought into the form 

Hb — £ uj n a n a n + const. (32) 

n#0 

by solving an eigenvalue problem, and the a n (a n ) are also boson creation (annihilation) operators. As Hb in Eq. 
(30) is a sum of commuting terms Hs, n only the coupling of the four operators b n , b- n , b n , b_ n has to be considered 
at a time. Momentum conservation further reduces the 4x4 eigenvalue problem to a 2 x 2 problem. The eigenvalue 
problem is obtained by looking at the Heisenberg equations of motion for b n (t) = e lHBt b n e~ lHBt . This involves the 
commutators 



[b n ,H B ] \ f k n (v F +g 4 (k n )/{2n)) fc„ ff2 (fc„)/(27r) \f b n 
[bl n ,H B }J \ -k n g 2 {k n )/{2ir) -k n (v F + g 4 {k n ) / (2ir)) J \bl 

The matrix M^") on the rhs of Eq. (33) has the eigenvalues 



(33) 




ut = ±v F k n Ml + ^) -f^A) = ±w „ (34) 



2ttv f J \ 2irv F 

As is not symmetric one has to distinguish left and right eigenvectors. We denote the left eigenvector to w+ as 

(c„, — s n ) where c„ and s n can be chosen as real quantities as MA™' itself is real. The formal solution of the Heisenberg 
equations of motion ib n — [b n ,H] etc. reads 



b n (t) \ _ _ iM (n) t ( b n 
bUt) ) { bl 



(35) 



7 



If we multiply this equation form the left with (c„, —s n ) we can replace the matrix by its eigenvalue (c„, —s n )e 
(c„, —s n )e~ luJrit which shows that the operator 



-!M<"'i 



a n = c n b n - s„6L„ (36) 

has the simple time evolution a n (t) = e~ lWn a n . With a proper normalization of the left eigenvector the a n are there- 
fore the boson operators which " diagonalize" Hb to the form Eq. (32). The a n should obey the boson commutation 

i„ the left eigenvector (c n ,—s n ) has to be normalized by 



relation [a„,a4] = 1. As Eq. (36) yields [a„,a4] = c\ 



1. Elementary calculation of the eigenvector then yields 



1 + gi(k n )/{2-KV F ) 



V(l + gA{k n )/(2nv F )Y - {g 2 {k n )/{2nv F )f 



(37) 



In the noninteracting limit g v — > the weights s n vanish and the "Bogoljubov transformation" in Eq. (36) reduces to 



The dispersion relation ujn^ — vp\k n \ of the bosons in the noninteracting case (see I) is in the interacting 



case replaced by uj Ti 
from Eq. (34) as 



v c {kn)\k n \ where the fc-dependent "sound velocity" usually called the "charge velocity" follows 



v c (k n ) 




gi{k„) + ff2(fcn) 
2itvf 



9A{k n ) - 92{k n ) 



(38) 



For our original model g^lk) = g 2 {k) = v(k) it is given by v c (k n ) = Vp(l + v(k n ) / {ttvf)) 1 ^ 2 ■ The corresponding 
dispersion is shown in Fig. 3 for the two potentials v presented in Fig. 1. 



4.0 





FIG. 3. Dispersion of the bosonic modes in Eq. (34) for the two potentials shown in Fig. 1. For the potential which is 
constant up to k c the dispersion is discontinous while for the smooth potential the approach to the noninteracting limit for 
large k is gradual. 

For \k n \/k c <C 1 the boson dispersion is approximately linear uj n « w c |fc„| with v c = v c (0). From Eqs. (31) and (38) 
we obtain the relation 

v c = Vmvj (39) 

first pointed out by HaldaneH. For large momenta \ k n \ 3> k c the dispersion of the noninteracting case is approached. In 
the hgure we have used a purely repulsive interaction (v(k) > 0) which leads to v c (k) > vf- For attractive interaction 
there is an instability if —v(q) gets larger than TrvpB, leading to phase separation. In the following we concentrate on 
repulsive interactions. 

The interacting ground state |{0} Q , N + , N- > with the integers N a = N a — mo — 1 is determined by 
c* t i|{0} Q , N + , 7V_ >= 0. It can be constructed by a unitary transformation out of the noninteracting ground state 
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\{0} b ,N+,N. >= TT cl_ _ TT c[ + |Vac > . (40) 



Here we have allowed the number of right and left movers to be different. Our original model corresponds to N + = 
N- = rip. In the following we calculate expectation values in theinteracting ground state by algebraic methods using 
an|{0} a , N+, N- >= 0. Therefore the explicit form of this sfcateQ is not needed. The excited states for given particle 
numbers N + , N_ are constructed as for harmonic oscillatorsQ 

n(rJ) mt 
Ud=-\{0} a ,N +: N_ > . (41) 

Their energy is given by E (N + ,N-) + J2i^to m i^i- 

The nature of the interacting ground state shows up clearly in the momentum distribution < n/ £n >. In the 
following we calculate < hk n , > and show how the Fermi step of the noninteracting limit is modified near the right 
Fermi point fop. As in I it can be obtained by first calculating < ^ + (u)ip + (v) >, where 4>+(v) is the auxiliary field 
operator defined in (B.6) and then using the inversion formula (B.24) to calculate the momentum distribution. 

With the use of the bosonization formula (B.25) the calculation is quite simple. As we want to use 
Q!n|{0}aj A+, A_ >= we express in the operator function 4> + defined in (B.9) the b n in terms of the a n ,ct_ n . 
From a n = c n b n — s n b^_ n in Eq. (36) one obtains the inversion 

b n = c n a n + s n o?_ n . (42) 

This yields 

i{<j)+{u) - <t>+{v)) = V — (c n a n + s n a'_ n ). (43) 

n—l v 

As the a n and a~t n commute, cxp[— i(<fi + (u) — (j>+(v))] can be written as a product of two exponentials with the 
annihilation part to the right. If applied to the interacting ground state this factor yields unity, i.e. 



-e:„( 



vt 



e -iW+(«)-^+("))|{ } a) iV + ,iV_ >=e ^"^V ^ J " _ "|{0} tt ,A>+,A>_ > . (44) 

The expectation value of V>+ («)■(/>+ (w) as seen from Eq. (B.25) involves the overlap of the above state with itself. 
Using the Baker Hausdorff formula e A e B — e B e A eS A,B ^ we obtain 

e -iN+(u-v) ( ^2s 2 ] 

< $.(uty+(») >= 1 _ ei(M _„ +i0) exp |- }2 ~f 0- ~ cos ^ u ~ ■ ( 45 ) 

This is exactly of the form as the finite temperature canonical expectation value for noninteracting fermions (Eq.(50) 
in I) but with the Bose function b(nAf3) replaced by which depends on the interaction as shown in Eq.(37). For 
arbitrary potentials therefore < n mj + > can be calculated by the recursive method outlined in the appendix of I. The 
sum in the exponent on the rhs of Eq. (45) is especially simple, if s 2 n decays exponentially 

The corresponding potential v(k n ) can be obtained using Eq. (37) by solving a quadratic equation. The factor 2 in 
the exponent on the rhs of Eq. (46) was introduced because v{k n ) is proportional to s n for large momenta. It is 
shown as the full line in Fig. 1 for s(0) 2 — 0.125. If the cosine is written as a sum of two exponentials sums of the 
type YlnLi Q n / n = ~ l n (l — l) appear in the exponent on the rhs of Eq. (45). This yields 



-iN+ (u — v) 



(1 - e- 2 ' n ' 



v+2i/n c ) \{\ — ^—i{u—v—2i/n c )\ 



(47) 



The calculation of < n„ li+ > simplifies if we go to the thermodynamic limit L —> oo, N+/L = const. Then it is 
appropriate to switch to the "physical" field operators ip a (x) defined as in Eq. (2) but with c n replaced by c„ iQ . They 
are related to the auxiliary field operators defined in (B.6) by 







The momentum distribution < hk m ,+ > follows from the inversion formula (B.22) and the fact that < ^ + (u)^)+{v) > 
is a function of u — v only as 

,L/2 

< h kn .+ > = / dxe ik " x < ipl(x)ip + {0) > . (49) 

J -L/2 

If one now takes the limit L — > oo, 2nN+/ L = kF one obtains for fixed x from Eq. (47) 

<^>"5^W'"' (50) 

In contrast to the noninteracting case this zero temperature time- independent " Green's function" for large x does not 
decay as (l/x) d , where d= 1 is the spatial dimension, but as (l/x) 1+2s i - \ Therefore the quantity 

a = 2s 2 (0) (51) 

is called the "anomalous dimension". The momentum distribution in the thermodynamic limit will be denoted n a {k). 
The behaviour of n+(k) near the right Fermi point is most easily studied by looking at its derivative 



dn+{k) 1 f°° ( {2/k c 



d~k 2n J_ QC] \x 2 + (2/k, 



* 2 (0) 



Ak x 



dx, (52) 



where k = k—kp- In the noninteracting case s 2 (0) = and (dn+(k) / 'dk)^ = —S(k), where the negative delta function 
is due to the step at k = kp- For finite values of s 2 (0) this singularity for k = is weakened. For 2s 2 (0) > 1 the 
function of a: Fourier transformed in Eq. (52) decays faster than l/\x\ and (dn+/dk)\z_ is finite. For < 2s 2 (0) < 1 
one obtains the behaviour mentioned in section II. As dn + /dk is an even function of k, it is sufficient to consider 
k > 0. With the substitution u — kx one obtains 



dn + mul w-i' 1 



2 (0) 



{2k/k c ) a - L ( )/ = cosudu. (53) 

dk ' ' ^ TtVo \u 2 + (2k/k c ) 2 J ^ ' 

For < 2s 2 (0) < 1 the limit k — > in the integrand can be performed and the integral is finite, but dn+/dk diverges 
at k = because of the prefactor of the integral. Integration with respect to k finally yields Eq. (16). At finite 
temperatures the derivative of n+ at kF diverges like for T — > 0. The special case a = 1 requires a separate 

treatment not given here. 

The power law behaviour Eq. (16) of n + {k) has been derived here for the special potential shown in Fig. 1 leading 
to the exponential decay Eq. (45) of s 2 . Different potentials with a different decay of s 2 (k n ) lead for \k\/k c <C 1 to 
the same qualitative behaviour as for our special potential if v(k) goes to constant value for k — > 0. This implies a 
constant value s 2 (0) and the anomalous dimension is given by Eq. (51). The anomalous dimension can be expressed 
in terms of the velocities vn and vj using Eqs. (31) and (37) 

2, 2 (0) = ^±^-l.^ + l-2), (54) 

y/VNVj 2 K 



where K = y/vj/vjy > is less than one for repulsive interactions and larger than one for attractive interactions. 

The power law behaviour of n+(k) for interactions v(k) with a finite limit v(0) in one spatial dimension is different 
from the " Fermi liquid" type behaviour of a usual three dimensional conductors in which the momentum distribution 
of the interacting ground states has a finite discontinuity of size < Z < 1 at the Fermi momentum! 2 ]. We note that 
for a finite one-dimensional system and small values of a the distinction between the power-law behaviour and a finite 
discontinuity is not straightforward. In Fig. 4 we show the momentum distribution for a step potential (see Fig. 1) 
with v(0)/(2itvf) — 1-5 and 0.3 corresponding to a = 0.25 and a = 0.0277 for different system sizes. 
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FIG. 4. Momentum distribution for the step potential of Fig.l with v(0)/(ttvf) = 3 and 0.6 corresponding to values of the 
anomalous dimension a — 0.25 and 0.0554. For the larger interaction the distribution is shown for different system sizes (stars: 
n c = 100, circles: n c = 200). For the smaller interaction the exact result for n c = 20 (squares) is compared to the result from 
perturbation theory (triangles). 

The size of the system is specified by the dimensionless ratio n c — (k c L) /2tt = (L/R)/2ir where R is the spatial range 
of the interaction. For the smaller interaction we also compare the exact result with the lowest order perturbation 
theory result Eq. (15) in the range < k — kp < k c . As the dimensionless coupling constant v(0)/(2ttvf) = 0.3 is not 
really much smaller than one the deviations are already well visible. 

Power laws occur for the TL- model also in other physical quantities. In the next section we present a .short discussion 
of one-particle spectral functions which largely determine the current in a photoemission experimentcnij . 



IV. SPECTRAL PROPERTIES AND PHOTOEMISSION 



For noninteracting electrons the "density of one- particle states per unit volume" (DOS) po(e) plays an important 
role for the calculation of thermodynamic properties!!!. In one dimension it is defined as 



Po(eJ 



(55) 



k„ 



In the thermodynamic limit X)fc(') ~~ * f(')dk and for a system of noninteracting right movers (e& = vpk) one 
obtains po( e ) = I/C^ttUf)) i.e. an energy independent DOS. For photoemission only the occupied states matter and 
one defines the "occupied DOS" Pq(c) by multiplying the delta function in Eq. (55) with the Fermi function /(efe„), 
which at T = is the step function at the Fermi energy tp. This simple DOS concept has to be generalized properly 
when the fermions interact. For a translational invariant system as discussed in the previous sections one defines 



p K (e) 



IE 



< MN)\cl5(e + H- E (N))c n \MN) > 



(56) 



as the generalization of the occupied DOS relevant for photoemissiorS. In the noninteracting limit c n \(f>o(N) > is the 
normalized Slater determinant where the state with momentum k n has been removed relative to the N-electron ground 
state. As its energy is given by E (N) — ek n the spectral function p < (e) reduces to the sum of delta functions as in 
Eq. (55) restricted to the occupied momenta. The individual terms in the sum of the rhs of Eq. (56) enter "angular 
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resolved photoemission"llSll3, while p < (e) is needed for the description of angular integrated photoemissionEl. For 
translational invariant systems expectation values with cjj replaced by c n , with n' ^ n in Eq. (56) vanish because of 
momentum conservation and one can write p < (e) in the form 

p<(e) = < MNM f (x = 0) (J- J e <(«+H-^W)*dt^ ^ x = 0)\cj> Q (N) > (57) 



= J <MNW(x = 0Mx = 0,t)\MN) >e ^ = I & < (t)e let — , 

where ^/>(x, t) = e t ip{x)e~ i is the field operator in the Heisenberg representation. If we switch from the original 
nonrelativistic model to the TL-model the spectral function acquires an additional a label. In the following we sketch 
the calculation of p$(e). For that purpose it is convenient to express the field operator ip+(v) Eq. (B.8) in terms 
of the a n , o4 instead of the b n , b^ n because the former have a simple time dependence. In b m = c m a m + s m a_ m the 

operators a m and a! m commute. Therefore e l ^+^ (and e 1 ^^) in Eq. (B.8) can be written as a product of two 
exponentials with the annihilation operators to the right. After using the Baker-Hausdorff formula once in order to 
complete the process of "normal ordering" (bringing all annihilation operators to the right of the creation operators) 
one obtains 

^ +{x) = (j^j e *l<*>e*+<»> (58) 

1 

with ix+(x) = V —j= 



jm 

m—l 



Eoo 2 i 



This is one of the most important formulas for the physics of one-dimensional interacting fermions. The time 
dependent operator ip + (x,t) follows from Eq. (58) by replacing a m and a_ TO by a m e~ luJmt and a_ m e~ Mmt and U+ in 
+ (see Eq. (BJ23^Lby^U + (t). Various kinds of time dependent correlation functions can quite simply be calculated 
using this resultOfl'ESTli. Here we restrict ourselves to G+(i) in order to obtain p+(e). 

As U+ commutes with the bosonic operators we first show that the particle number changing operators U\_ and U+(i) 
lead to a simple time dependent factor. Using (B.17) one obtains 

Ule im U+e- im \{0} a ,N+,N- >= e^^+^-^+- 1 '^ t \{0} a ,N+,N- > . (59) 

As ip+(x) in Eq.(58) is normal ordered in the a's one has to use the Baker-Hausdorff formula only once to normal 
order ipl_(x = 0)ip+(x = 0,t). This yields 

ie^G<(t) = ^!M e [x + (x=o,t=o), x ;(x=o,t)] (60) 
L 




where p, = Eq(N +) iVL) — Eq(N + — 1,N~) is the chemical potential, which we "eliminate" by discussing p^_ as a 
function of e = e — p. The spectral function can be calculated analytically in the "universal" low energy regime 
— fc c min (vf,v c ) < e < for the step potential which has a strictly linear boson dispersion u n = v c \k n \ for \k n \ < k c 
(see Fig. 2) and a different linear dispersion w„ = t>F|fcn| for \k n \ > k c . This implies that we can write G+{t) as a 
product of two power series in exp(i2Trvpt/ L) and exp(i2irv c t/ L) 

iLe^G^t) = A 2 (n c ) ll + a£< W^*) ) (|^) e M^)). (61) 

\ m=n c + l / \n=0 / 

(n ) ~ (n) 

Now the Fourier integral can be performed analytically. For VFk c < e < the % and for u c fc c < e < the d m with 
?7i > n c do not contribute and one obtains in the universal regime — k c mm(vp, v c ) < e < 
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Lp<(e) = A 2 (n c ) ■ £ d^S (i + n^v c ) . (62) 



n>0 



The d$h"^ for n < n c can be given analytically as the sum in the exponent in Eq. (60) is of the form (1 + a) ■ X)"=i 2 ™/ n 
with z = exp(i2irv c t/ L) and a = 2s 2 (0). For the calculation of the Sh"^ for n < n c one makes no error to extend 
the sum to infinity, which yields exp [(1 + a) ^2 n z n /n] = (1 — z)~( 1+a \ Now one can use the well known power 
expansion of this function. Alternatively one can show that the recursion relation in (A4X of I can simply be solved 
explicitly. If we also use X]™=i V n ~~ * ui7i c + C, where C = 0.57721 ... is Euler's constantly we find the exact weights 
of the delta peaks in Eq. (62) 

— aC n 

for n c > n > 1 and cIq 1 ^ — 1. For n c > n 1 one can use the product representation of Euler's Gamma functionEl 
to write the weights as 

A 2 (n c )d^= T^rr-T (-V(l + 0(l/n c )) . (64) 
r(l + a) \n c J 

Now the thermodynamic limit can be performed in Eq. (62) and one obtains with e c = v c k c for e < p, in the universal 
regime 

The spectral weight near (below) the chemical potential \l is suppressed by the same power law behaviour as in 
n + (k) — 1/2. This type of power law occurs also for potentials with a finite v(0) but a different decay for momenta 
large compared to k c . But then the prefactor is different and the pojwer law behaviour does not hold exactly in a finite 
range as Eq. (65) but only asymptotically in the limit (— e)/e c — > OlZI. This power law suppression of spectral weight 
was claimed to be expjerimentally verified in photoemission experiments with an organic compound which consists of 
weakly coupled chainalS. 

Reliable results for more complicated models of interacting fermions can often only be obtained using numerical 
methods for systems of rather small size. As our result Eq. (63) is exact also for finite systems we shortly discuss 
how one can extract the anomalous dimension even for small values of n c . As Eq. (64) only holds for n 1 it is not 
useful to fit the weights to the power law form. It is much easier to look at the scaling behaviour of the peak at e = [i 
with weight A 2 {n c )d!^" : ^ — e~~ aC / (n c ) a with system size. From a log-log-plot one can directly infer the value of a. 

For — e ^> e c the spectral function p+(e) approaches the poninteracting value 1/(2tvvf). The crossover from the 
power law behaviour to this high energy limit is discussed ir£j. 

Finally we want to point out that the low energy power law behaviour of /0+(e) given by Eq. (65) holds for 
arbitrary values of a. No additional linear contribution occurs for a > 1 as in the case of the momentum distribution 
n + (k) — 1/2. This difference is easily lost without a proper treatment of the cut-offs. 



V. SPIN-CHARGE SEPARATION IN THE MODEL INCLUDING SPIN 



Electrons are spin one-half particles and for their description it is necessary to include the spin degree of freedom 
in the model. We choose a fixed quantization axis and denote the two spin states by a —] 1 [. Then the field operator 
tp(x) acquires a spin index as the c n . The density operator of spin projection a is defined as p a {x) = ipl.(x)ipa(x) and 
in the two-body interaction in Eq. (6) V(x — x')p(x)p(x) is replaced by J2 a a 1 Varf ( x ~ x ')p<t( x )p<t' {x'). Usually the 
interaction V aa > (x — x') is independent of the spin indicesi. After the linearization and introduction of the right and 
left movers the p n ^ a in Eq. (20) obtain an additional spin index, as well as the boson operators in Eq. (27). The 
TL-Hamiltonian for spin 1/2 then reads 



H, 



1/2) 



TL 



E E 



k n V F S a 
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fcnfflf 7 (fcn) / ii 

2^ S 



I 



n,cr — n,a 



9%" (0) 



, +b. 



n.ij' u na 



(66) 



2tt 



9T 



(0) 



The interaction matrix elements g^ 17 (k) are taken as depending on the relative spin only, i.e. g1° and g°' a are 
independent of the spin label a. The interaction term in Eq. (66) couples the electrons with different spin. With the 
assumption on the g° " it is possible to switch to new boson operators b n ^ a with a = c, s 



K, c = -/= (K,i + K,i) 



(67) 



which obey [b a 

(c) and "spin" 
elements g v>a (.q) via 



<^aa' ^nn'l. The kinetic energy can also be expressed in terms of "charge" 



(s) Bose operators using b\ j&n.T 



+ &„|&ni = bl hC b n , c + b^ gbn^. If we define the interaction matrix 



and define AC 



c(s) 



9",c{q) = gl a {q) + 9° 
9vAi) = 9™{q) -gl 

(W T ±Ni)/V2 we can write H^[ 2) as 

r(i/2) 



'(q) 



H 



TL 



Htl.c + Htl„ 



(68) 



(69) 



where the Htl,o, are of the form Eq. (30,31) but the interaction matrix elements have the additional label a. The 
two terms on the rhs of Eq. (69) commute, i.e. the "charge" and "spin" excitation are completely independent. This 
is usually called "spin-charge separation"B. The "diagonalization" of the two separate parts proceeds exactly as in 
section III and the low energy excitations are "massless bosons" ii)„ j0 = i> a |& n | with the charge velocity v c and the 
spin velocity v s usually different. For a spin independent interaction g„ jS in Eq. (68) vanishes and v s equals the 
Fermi velocity vf- As discussed in the next section the Hamiltonian H^/ 2 ^ with properly chosen interaction matrix 
elements describes correctly the low energy physics of more complicated Id models e.g. on a lattice. Then v s and Vf 
are generally different even if the interaction on the lattice is spin independent. 

The bosonization of the field operator ip atiy (x) proceeds exactly as in appendix B. In order to calculate quantities 
like the momentum distribution {fik.+^) it is then useful to express the 6„ j(T in i<j) +ja -(v) like in (B.9) in terms of the 
charge and spin bosons 



1 

71 



E 



-{br, 



'v) ± ic 



,(«) 



(70) 



where the plus (minus) sign holds for the up (down) spin. Therefore e.g. the function (tp' + |(a;)^+ 1 -|-(0)) is a product 
of a charge and a spin part. The separate factors are calculated as in the spinless case. The only difference is a factor 
1/2 in the exponent due to the additional factor l/y/2 in Eq. (70). The final result is exactly of the form Eq. (45) 



but with 2s„ replaced by s 



2 

n.c 



This leads to the power law behaviour presented in Eq. (16) as in the spinless 



case but the anomalous dimension given by 



a = s 2 c (0) + s 2 s (0) = a c + 



(71) 



,1/2 



The individual contributions can be expressed in terms of the K a = {vj. a /vN, a 

(t) needed in the calculation of the integrated spectral function p+ a (e) takes the 



to Eq. (54). Also the result for G< 
form of a product of a "spin" and a "charge" factor. In each factor the weight 1 



as a a = (K a - l) 2 /{AK a ) similar 

< 



2s„ in the exponent of Eq. (60) is 



replaced by 1/2 



Therefore the Fourier transform of the individual factors are in the low energy regime of the 



power law form Eq. (65) with a replaced by — ^ + a a . The convolution of the two power laws finally for e < /i yields 
CT (e) ~ (jU — e) Q in the asymptotic regime as in the spinless model. Alternatively this follows more directly from 
the fact that Gf a (t) decays as t~ t - 1+a > for large times. 

In order to see a drastic difference between the model with and without spin, one has to calculate the individual 
terms in the sum in Eq. (56). The delta-^pak-Si of the noninteracting model are broadened into one "infrarsd," r-ppwer 
law threshold in the model without spinOJI13'til and two power law singularities in the model including spinli3ii-l 
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VI. THE LUTTINGER LIQUID PHENOMENOLOGY 



We have seen that the low energy physics of the TL- model is determined by two parameters v c and K or vn and vj 
in the spinless model and four parameters (v c , K c ; v s ,K s ) or (wjv c , vj c ,vn s , vj s ) in the model including spin. They also 
determine the low temperature thermodynamics of the model. We shortly discuss the model including spin. The low 
temperature specific heat of a linear chain is gj«en by 7T, where the constant 7 is inversely proportional to the sound 
velocity. This is the Id version of Debye's lawu'Q. The ratio of 7 and the corresponding value 70 for the noninteracting 
fermions is therefore 

^ -\ £ + v, ) ■ 

At zero temperature we can express the pressure p as p — — (OEq/DL) n and the compressibility k = —{DL/dp)iq/L = 
(L{d 2 E /dL 2 ) N y 1 = ((N 2 /L){d 2 E /dN 2 ) L y\ Because of Af c = (A/| +J\f l ) /y/2 the TV-dependence of the ground 
state energy of the TL-model using Eq. (30) is given by v/v c (7r/4L)iV 2 leading to 

/ d 2 E \ 7T n v c 



L \jN^)r 2 VN ^ 2 1£ (?3) 



The compressibility ratio is therefore given by 



k/k q = — K c . (74) 



The factor 1/2 in Eq. (73) is absent in the spinless case. We leave it as an exercise to show that one obtains similarly 
for the spin susceptibility 

X/Xo = —-K s . (75) 
v s 

These relations show that the K a which determine the contributions a a to the anomalous dimension a and therefore 
the power law behaviour of the Green's functions also enter the low temperature thermodynamics. 

It was an important insight by Haldanetl to realize that the low energy physics, of the TL-model is generic for 
interacting fermions in one dimension. In the language of the renormalization groupta the TL-Hamiltonian is the fixed 
point Hamiltonian for Id fermions with a repulsive interaction. This is the main reason for the importance of the 
TL-model for the physics in one dimension. 

Usually the LL parameters K a ,v a of the effective TL-model which describes the low energy physics of a given 
model can only be calculated approximately. An example is the original Tomonaga model with the quadratic energy 
dispersion (with or without spin) where it is not assumed that the interaction cutoff is much smaller than the 
Fermi momentum, i.e. an arbitrary repulsive interaction is used. An exception is the extreme short range limit 
V(x) — Vo8(x). For the spinless model V then vanishes because the Pauli principle forbids two fermions at the same 
space point (?p(x)?p(x) — in Eq. (1)) . For the corresponding model including spin ground state properties can be 
calculated exactly using the Bethe-ansatz methodH. The ratios K c /v c and K s /v s of the effective TL-model which 
describes the low energy physics of the model can be inferred by calculating the compressibility and spin susceptibility 

using Eqs. (74) and (75). The charge and spin velocities can be obtained as {E c [s) {2n / L) - E )/(2tt/L), where 
E^ s \2tt/L) is the lowest energy for charge (spin) excited states with total momentum 2n/L. From this information 
one can obtain the anomalous dimension, which determines the decay of the Green's function without the need to 
calculate this function. For Id models on a lattice, which cannot be solved exactly by the Bethe-ansatz method one 
can determine the LL parameters approximately by numerically calculating the quantities discussed above for small 
systems. In the case of weak interactions perturbation theory can be used. We show this to leading order for the 
spinless Tomonaga model with a short range interaction. The ground state energy up to first order in the interaction 
is given by the expectation value of H in the Fermi sea \F > 



E {1) 

/ - 2m e y L ' 2L 



,L t - ^ ~ lip Tip 

E ^r(-?) 2 + W E E (v(0)-v(k n -k' n )). (76) 



The interaction contributes the Hartree and the Fock term. In the large system limit the sums can be replaced by 
integrals in the standard way and one obtains 
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( r) 2 F {1) \ 

L [m^J =™F+v{Q)-v(2k F ) = K(v c lK)^. (77) 

In order to obtain the charge velocity to linear order in the interaction we have to calculate the expectation value of 
H in the state |^ >= c^ jr+1 c„ j; , \F >. Taking the difference with Eq and dividing by 2tt/L yields 

(« c ) (1) = v F + (v(0) - v(2k F ))/(2n). (78) 

The contribution of the interaction is due to the different Fock energies in the states \=r- > and \F >. Combining the 
results from Eqs. (77) and (78) gives K = 1 — (v(0) — v(2k F )) / (2irv F ) + ... which leads to 

(2) 1 f v(0)~v(2k F ) 

a = - \ 

2 V 2ttv f 

If one generalizes the calculation of < h m >' 2 - ) in section II to the case of a short range interaction one obtains the same 
result for the anomalous dimension. This shows to leading order perturbation theory that the low energy physcis of 
the Tomonaga model with an arbitrary range interaction is that of the TL model with 32(0) = 34(0) = w(0) — v(2k F ). 

Because of the limited length of the paper the rapidly expanding field of electronic transport in one-dimensional 
interacting Fermi systems was not discussed. Important insights can rather quickly be obtained from the simple 
exercise-jOf normal ordering of ip\_(x)ip-(x) with the use of Eq.(58) and the corresponding expression for the left 
moveraij. The interested reader should look up the relevant references in the review article citedo, especially the 
important work by Kane and FisherE2l. 
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APPENDIX A: THE KRONIG IDENTITY 



In this appendix we generalize an identity found by Kronig in the context of the (now long forgotten) "neutrino 
theory of light"liJ to the case of a finite cut-off ks- For the right movers it reads 

oo oo _^ 

Y lh \ bl= Y (n* + *no)<4,+Cm,+ - 2 • (A- 1 ) 

1=1 m= — m 

For the proof we simply evaluate the lhs using Eqs. (20) and (27) 



oo oo 



Y l tf bi = Y Y c^ + ; i+ c n , + 4 i+ c n+i , + + y Y 4+i,+ c «.+ c t,+ c n'+i.+ 

1 = 1 1=1 n= — m 1 = 1 n^n' 

= ii+i 2 , (A.2) 

where we have split the double sum over n and n' in the contributions from n = n' and n n' . We will show that 
the operator A\ yields the rhs of Eq. (A.l) while A^ vanishes 

00 00 00 00 



Ai=Y Y ~ Y Y (A.3) 

1 = 1 n= — m 1=1 n= — m 

= Y ( m + m o)n m ,+ - 2^+~ ■ A/ +) ' 

m--mo 

The second equality follows from a change of summation variables and using , = n m ,+ - For the discussion of A2 
we split the double sum into the contributions from n > n' and n < n' 

^2= y \Y c i+i, + c ",+ c i',+ c «'+i,+ + Y c i+i,+ c ™.+ c i',+ c ™'+'-+ [ ■ ( A - 4 ) 

1=1 \n>n' n<n' ) 

In the first term we write n = n' + I' and use /' > 1 as the new summation variable and in the second term we write 
n' = n + V. This yields 

oo 00 ( oc 00 ~\ 

^2 = Y Y \ Y C l+l+l', + C l', + C n'+l, + Cn'+l\+ + Y C i+i, + C i+i', + C «+'+'', + C n,+ > ■ ( A ' 5 ) 
1=1 l' = l Kn' = — m n= — m ) 

As the products of the four operators are both antisymmetric in I and I' the operator A2 vanishes after performing 
the summations, which completes the proof. If we multiply Eq. (A.l) with vf{2tt/L) and use the analogous relation 
to (A.l) for the left movers we obtain Eq. (29) for the kinetic energy Ttl- We note that in the proof no use is 
made of the commutation relations of the 6; and b\ . The Kronig identity holds as an exact operator identity without 
performing the limit ks — *■ 00. 



APPENDIX B: BOSONIZATION OF THE FIELD OPERATOR 



In this appendix we generalize chapter V. of I on the bosonization of fermionic operators to the case of a single 
fermion operator . A- new straightforward construction of the particle number changing part is given which simplifies 
earlier presentations^. 

We start with a simple relation which holds for operators 6, W obeying [b, b'] = 1. It reads 

[6,e Abt ] = Ae Ab+ . (B.l) 

To prove it one differentiates 6(A) = e~ Abf 6e Afct with respect to A and obtains db(X)/dX = 1. Therefore 6(A) = 6 + A1 = 
e~ Ab 6e Afct . Multiplying from the left with e Abf yields (B.l). Next we define the following operators linear in the boson 
operators b n , 6^ 
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A + = Y,^nbl ; B_ = Y,»nbn (B.2) 

with arbitrary constants A n and [i n . Using (B.l) and elementary generalizations one can prove the commutation 
relations 

[b m ,e B -e A +] = \ m e B -e A + (B.3) 
[bl,e B -e A +] = -v m e B -e A + . 

In the following we consider fermionic operators S which obey [b m , S] = —X m S and [6^, 5] = n m S. Then the operator 
O = Se B ~e A + commutes with all b m and b^. We therefore write 

S = Oe- A +e- B -, (B.4) 

and subsequently construct O such that both sides of Eq. (B.4) yield identical matrix elements. 

In the following we present the bosonization of the field operator for the right movers described by the c; i+ and just 
mention the corresponding result for the left movers. The C{,+ obey for to > the following commutation relations 

[b m ,ci.+\ = ]=c l+rn . + , [b f m ,c L+ } = — ci- mt+ , (B.5) 

V m v m 

The auxiliary field operator ip+(v) defined as 

oo 

$+{v)= lim Ve^ Ci .+ = V e av d i+ . (B.6) 

l — — oo 

obeys commutation relations as the operator S discussed above 

[b m ,j>+(v)] = -^=e~ im ^ + {v) ; [blMv)} = -- L e ™> + («) . (B.7) 



/TO 

Therefore it is of the form 

$+(v) = 6 + {v)e l ^+ {v) e l ^ {v \ (B.8) 

where the operator i<f>+{v) is given by 



n—1 * 



(B.9) 



As ip+(v) reduces the number of right movers by one, the operator 0+(v), which commutes with all b m , b^, also must 
have this property. In order to determine 0+(v) we work with the eigenstates of the noninteracting system (compare 
Eq. (41)) 

\{ mi } b ,N + ,N_>=]] y -jL=-\{0} b ,N + ,N_> . (B.10) 

After the limit m — > oo the N a can run from minus to plus infinity. It is easy to see that O + (v)\{0} b} N + , N- > has 
no overlap to excited states 

< {mi} b ,N+ - l,N-\6 + (v)\{0} b ,N+,N- > = 

<{0} b ,N + -l,N_\]J y -jL=d + (v)\{0} b ,N + ,N- > . (B.ll) 



i 



'mil 



As 0(v) commutes with the bi the rhs of Eq. (B.ll) vanishes unless all to/ are zero. This implies 
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O + {v)\{Q} b ,N + ,N->=c+ + j_(v)\{0} b ,N + -l,N-> , (B.12) 

where ct . (v) is a c-number. In order to determine the . we calculate ({0}b, N+ — 1, A r _|-0(w)|{O}h, 7V + , iV_) 

using (B.6) as well as (B.8). In the calculation of the matrix element with the fermionic form (B.6) we use Eq. (40) 
which yields 

({0} b ,N+-l,N.\c l , + \{0} b ,N + ,N_) = (-lf-d l ^ + . (B.13) 

The factor (— 1) occurs because we have to commute q.+ through the product of 7V_ fermionic operators of the 
left movers if we assume mo to be odd. We note that no such factor occurs for the corresponding matrix element of 
the left movers. The calculation of the ground state to ground state matrix element of tp+( v ) using (B.8) is simple as 
both exponentials involving the boson operators can be replaced by the unit operator and the matrix element is just 
c~~ . The comparison therefore yields 

4 +> ff>) = (-l)*-e to * + (B.14) 
and ^ (v) — e~"' N - . Together with (B.12) and the definition N a = M a — {mo + 1)1 this implies 



O + (v)e- l "+ l >(-l)"-\{0} b ,N + ,N_) =\{0} b ,N + -l,N.) . (B.15) 

If we apply the operator 0+(u)e~' iA, + l '(— 1)-^- to the states in Eq. (B.10) and use again that 0+(v) commutes with 
the boson operators we obtain 

d + (v)e- i *+ v (-l)*-\{m l } b ,N + ,N^) = \{ mi } b , N+ - 1, N-) . (B.16) 

This shows that the operator U+ = d+(v)e~ lJ ^ +v (—l)^- is independent of v and given by 

U+= E \{mi}b,N+-l,N-)({m l } b ,N+,N-\ . (B.17) 

N+,N- i m i} 

It follows immediately that U+ is unitary, i.e. U+U\_ = U\U+ = 1. From (B.15) one can infer that for arbitrary 
functions / of the number operator N+ one has U+f(M+) — /(A+ + 1)E/+. Concerning the ./V + -dependence the 
operator U + resembles a one-dimensional tight-binding nearest neighbour hopping operatorH. Therefore the eigenstates 
of U + are just the corresponding "Blochstates" 

\{mi} b ,k,N-)= -= V e**+|{m,} 6 , (B.18) 

V 27T - 

JV + = -oo 

and the eigenvalue of U+ is given by e tk . If one defines the operator 

k+= [ \{mi}b,k,N-)k({ mi } b ,k,N-\dk (B.19) 

the operator U+ can be written as 

U+ = e ik + . (B.20) 

In the "^-representation" the particle number operator A/+ acts like the differential operator id/dk, which seems to 
imply the commutation relation [A/+, fc+] = it ■ This is not correct as an operator identity in the full Hilbert space. 
In the fc-representation the operator fc+ acts as a multiplication operator which destroys the 27r-periodicity in k unless 
the state fulfills ({mi} b ,k = tt,7V„|^) = = ({mi} b ,k = -ir,N-\ip). 
The direct calculation of the commutator using 

Af+= Yl E \{mh,N+,N-)N + ({mi} h ,N+,N-\ . (B.21) 

N+,N- {mi} 
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yields 



[tf + ,k+]=i 1- 2 ^H J^\{mi} b ,k = 7r,N_)({m l } b ,k = TT,N-\\ (B.22) 

\ {mi} JV_ / 

in contrast to statements in the literature!. 
To summarize we have shown that 

6+{v) = U+e^+^-l)*- = e^+e^+^-l)^"- (B.23) 

In the corresponding expression 0_(u) = e' lk - e ~ lJ ^- u no factor (— 1)-^+ appears and therefore 0+(v) and 0-(u) 
anticommute, which is necessary to enforce anticommutation relations between ip+{v) and tp-(u). The bosonization 
of q ; + itself is obtained by the inversion formula to (B.6) 

ci,+ = ^ / e- M i,+ (v)dv . (B.24) 



— 7T 



We finally show how the bosonization of a single fermion operator described in this appendix relates to the bosonization 
of the product ip+{u)''4 , + ( v ) described in I. Using (B.8) and (B.21) we can write the product as 

(w) = e -^lWe- l * +(u) e^+ {v) e l *+ iv) e- lji:f +( u - v) 

= e -i(0+(«)-0+("))g-i(0+(«)-0+(t')) e [0+(«),0 1 f (v)] e -iAf + (u-v) _ (B.25) 

Here we have used the Baker- Hausdorff formula e A e B = e B e A e^ A ' B \ valid if [A,B] commutes with A and B. The 
commutator in the exponent is given by 

-, 00 p inu p —in'v 00 

n,n / — 1 n— 1 

= - lim [in (l - e <(*-«+*?)^ . (B.26) 

In order to make the infinite series convergent we have added a convergence factor e~ r ' n with -q going to zero in the 
end of the calculation. This is necessary^s we had to perform the limit mo — > oo in Eq. (B.6), i.e. we work with an 
infinite Dirac see as in Luttinger's papera. Combining Eqs. (B.23) and (B.24) gives 

= e-*W+M-*+(«" e -^ + («)-<M»» (B.27) 

which corresponds to Eqs. (47) and (49) of I. 

A convergence factor also has to be introduced explicitely in the bosonization formula for the field operator (B.8) 
if no normal ordering is used. If one uses the Baker-Hausdorff formula and (incorrectly) assumes the commutation 
relation [A/+,fe+] = il one obtains 

&,(«) = lim — e * ] ("), (B.28) 

where the "phase field" operator 

§W(u) = k a + a(Af a - l)u + ^ ] (u) + (^(u)^ + |(1 + a)Af. a , (B.29) 

with 

p inu 

</>M( M ) = -iJ2®(an)-= b n e"!"!"/ 2 (B.30) 



n#0 

contains the boson as well as the particle number changing contribution. It is used in various applications! e.g the 
discussion of transport. It is obviously more rigorous to use of the normal ordered form (B.8), as always done in the 
present paper. 
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